Exact relativistic tritium /3-decay endpoint spectrum in a hadron model 
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We present the relativistic calculation of the /3-decay of tritium in a hadron model. The elementary 
particle treatment (EPT) of the transition :i H ^ 3 He + e~~ + V e is performed in analogy with the 
description of the /3-decay of neutron. The effects of higher order terms of hadron current and 
nuclear recoil are taken into account in this formalism. The relativistic Kurie function is derived 
and presented in a simple form suitable for the determination of neutrino masses from the shape of 
the endpoint spectrum. A connection with the commonly used Kurie function is established. 
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I. INTRODUCTION 



trino is a Majorana or a Dirac particle. 
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Neutrinos are one of the most intriguing and fascinat- 
ing fundamental particles, which make up the Universe. 
However, they are also one of the least understood par- 
ticles. Studies of neutrinos have played a crucial role in 
the understanding of elementary particle laws and their 
interactions. 

Three types of light neutrinos are known. The recent 
observation of neutrino oscillations P, 0, H, d, Q has now 
beyond doubt established the non-zero masses of neutri- 
nos, the flavor change and neutrino mixing. It has opened 
a new excited era in neutrino physics and represents a big 
step forward in our knowledge of neutrino properties and 
serves as solution of many problems in cosmology, ele- 
mentary particle physics, and astrophysics. 

While neutrino oscillation experiments are sensitive 
only to differences of squared neutrino masses, the 
neutrino mass measurements with tritium (Qp( H) = 
18.6 keV) and rhenium (Qp( lfi7 Re) = 2.47 keV) 0- 
decays yield direct information on the absolute neutrino 
mass scale. The idea underlying the measurement of neu- 
trino mass is actually fairly obvious. A long time ago, it 
was already pointed out by E. Fermi Q that the shape of 
the electron spectrum in nuclear /3-decay, near the kine- 
matical end point, is sensitive to the neutrino mass. 

Attempts to evaluate the rest mass of the neutrino ex- 
perimentally were already being undertaken long ago. In 
1940 one of the first kinematical measurements of neu- 
trino mass was performed by Hanna and Pontecorvo Q 
with a proportional chamber filled with tritium. A limit 
of ~ 1 keV on the neutrino mass was obtained, which 
was determined by the resolution of the detector. The 
Mainz Q and Troitsk Q tritium /3-decay experiments us- 
ing the magnetic adiabatic collimation technique, place 
the present upper limit on the mass of the electron neu- 
trino of 2.3 eV and 2.2 eV, respectively. The best pub- 
lished calorimetric limit to the electron neutrino mass 
obtained from the /3-spectrum of 187 Re is 15 eV [Toj ] . We 
note that the bounds on neutrino mass imposed by the 
shape of the spectrum are independent of whether neu- 



A next-generation tritium /3-decay experiment is the 
KArlsruhe TRItium Neutrino experiment (KATRIN) 
[l2l [l3l [l4T ] , which is presently in construction phase (It 
is planned to take data starting 2010). This experiment 
is projected for measurement of the neutrino mass with 
a sensitivity of 200 meV, which will have important im- 
plications for the theory of neutrino masses. If the result 
will be positive, it will imply a degenerate spectrum of 
neutrino masses. On the other hand, a negative result 
will be a very useful constraint. There is also a chance 
that the planned MARE experiment [ll[ based on ar- 
rays of rhenium low temperature microcalorimeters will 
be able to achieve sensitivity lower than 0.2 eV in fu- 
ture. The MARE approach would have totally different 
systematics with respect to the KATRIN. 



In view of an enormous experimental progress in the 
field there is a request for a highly accurate theoretical 
description of the electron energy spectrum in the de- 
termination of the neutrino masses from the shape of 
the endpoint spectrum. The subject of interest has been 
molecular effects in tritium beta decay [l5T | , radiative cor- 
rections [l6| . Lorentz invariance violations [l7T |. interac- 
tions beyond the standard model 1181 . relativistic form 
for the /3-decay endpoint spectrum [l9l[2f3| etc. 



The aim of this paper is to derive the relativistic form 
for the /3-endpoint spectrum in a hadron model. We shall 
take advantage of the fact that the nuclei 3 H and 3 He 
are, respectively, the nuclear analogs of the neutron and 
the proton, i.e., they form an isospin SU(2) doublet. A 
correspondence to the commonly used formulae will be 
established. We note that the considered approach is 
known also as Elementary Particle Treatment (EPT) of 
weak processes, which was developed by Kim and Pri- 
makoffEDl. 
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II. THE NUCLEAR PHYSICS DESCRIPTION 
OF TRITIUM /3-DECAY 

By neglecting neutrino mixing for simplicity and tak- 
ing into account only left-handed weak interaction, the 
electron energy spectrum for tritium /3-decay is 



N ( E e) = HT- 



dT 
d~E e 

GFV ^\M.E.\ 2 F(Z,E e ) Pe E e 



2tt 3 



x (E - E e ) y/(E - E e f 



(1) 



where Gp is the Fermi constant and V u d is the element 
of the Cabbibo-Kobayashi-Maskawa (CKM) matrix. p e , 
E e and Eq are the momentum, energy, and maximal end- 
point energy (in the case of zero neutrino mass) of the 
electron, respectively. F(Z, E) denotes the relativistic 
Coulomb factor. 

The transition is superallowed, a mix of Fermi and 
Gamow- Teller transitions. The absolute square of the 
nuclear matrix element is given by 



\M.E.\ 2 = ti\M F \ 2 + f 2 A \M GT f 



(2) 



where the Fermi and Gamow- Teller matrix elements take 
the form 



M F = < 3 He\J24\ 3 H>, 

fc=i 

3 

M GT = < 3 He^r+a^H > 



(3) 
(4) 



fc=i 



fv and j a are the vector and the axial-vector coupling 
constants of the nucleon, respectively. We note that the 
derivation of the differential decay rate in (|TJ) involves 
non-relativistic approximations and that only the S\ii 
states of outgoing leptons are taken into account. 

The Fermi matrix element can be evaluated by assum- 
ing the exact isospin symmetry as well as the fact that 
3 H and 3 He form an isospin doublet (T — 1/2) (the pro- 
jection T z = 1/2 is assigned to the 3 He and T z = —1/2 
to the 3 H) with the result Mp = 1. 

The absolute square of the Gamow- Teller matrix ele- 
ment can be deduced from the Ikeda sum rule by tak- 
ing into account that the Gamow- Teller operator has no 
radial dependence and thus can not scatter into higher 
shells. In 3 He the Is neutron level is already occupied 
by two neutrons and therefore in the transition p to n 
the neutron would need to be scattered into a higher or- 
bit (e. g., 2s) in the continuum, which is forbidden for 
the Gamow- Teller operator. Thus only 3 H -^ 3 He but 
not 3 H — > 3n can contribute to the Ikeda sum rule. In 
addition, there are no excited states of 3 He. As a conse- 
quence \Mqt\ 2 = 3. This result is in a good agreement 
with the recommended value \M GT \ = \/3(0.962 ± 0.002) 
obtained in nuclear structure calculation 241. 



The conserved vector current (CVC) hypothesis pro- 
posed by Feynman and Gell-Mann suggests that the vec- 
tor coupling constant fv is not renormalized in the nu- 
clear medium, i.e., fv — 1-0. The accurately measured [5- 
decaylifetime of tritium (T 1/2 ( 3 H) = 12.32±0.03 years) 
[2211231] is used to adjust the value of axial- vector coupling 
constant f a via the calculation of the theoretical half-life 



rp \-l/2 
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11T2 



Eo—rriu 



N{E e )dE e 



(5) 



In the computation of the integral over the electron en- 
ergy E e we adopted the relativistic Coulombic factor 
F(Z,E) [2^], which take into account the finite size of 
the nucleus. For \M GT \ 2 = 3 we found \f A \ = 1.247. The 
very good agreement between this result and the bare 
nucleon value \f A /fv\ = 1-2695 ± 0.0029 [H[ suggests 
that the axial-vector coupling constant is only weakly 
quenched in the tritium. 

The dependence of spectrum shape on the mass of neu- 
trino m„ in |1| follows from the phase volume factors 
only. The traditional way to look at the /3-spectrum data 
is to make a Kurie plot, where 



K(E e 



N(E e 



V F(Z, E e ) Pe E, 
GfV u< i 




(6) 



For zero mass neutrino, if K(E e ) is plotted against E e , 
the result is a straight line that crosses the E e axis at 
E e = E . For m u ^ the endpoint shifts to E max = 
Eq — m v and the rate near the endpoint is depressed, 
namely the Kurie plot has a kink at the endpoint. This 
distortion will be washed out at the experiment unless 
the energy resolution is comparable to m v . 

There are open questions related to the presented con- 
ventional approach for kinematical study of the /3-decay 
endpoint of 3 H . In particular, it is not known what the 
consequences of the considered non-relativistic approxi- 
mations are. Further, the effect of the nuclear recoil is 
not taken into account. It is also worth mentioning that 
the relativistic expression for the maximal electron en- 
ergy 



E" 



1 



2M f 



[M 2 +m 2 e - (M f + m v f 



gives a value about 3.4 eV lower than the considered 
approximation E 7 e nax ~ Mj - M f - m v [20] (M h M f 
and m e are masses of the tritium atom, 3 He + and the 
electron, respectively) . In view of the planned sensitivity 
of ~ 0.2 eV of the KATRIN experiment, there is a request 
for a consistent relativistic description of the /3-decay of 
tritium [2(J. 
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III. RELATIVISTIC /3-DECAY KINEMATICS IN 
HADRON MODEL 



We shall study the /3-decay of tritium, 

3 H -> 3 He + e~ +V e , (8) 
in an analogy with the /3-decay of a free neutron, 



P 



(9) 



as the spin-isospin characteristics of 3 H ( 3 He) nucleus 
and neutron (proton) are the same. The kinematics of 
the two processes above differ mostly due to different de- 
values and the Coulomb corrections. 

The invariant /3-decay amplitude is given by 

M = u(PMl- 75 )v(P u ) 



V2 
xu(P f ) 



G M {q 2 



G v (q 2 h a + i 



G A (q'h a l 5 ~G P (q')q a l5 u(P). 



(10) 



Here, q a = (Pf — Pi) a = (P e + P v )a is the momentum 
transferred to the hadron vertex. Pi — (Mi,0), Pf — 
(Mf, p f ), P e = (m e ,p e ) and P v = (m„,p„) are four 
momenta of the 3 H, 3 He, electron and antineutrino in 
the laboratory frame, respectively. 

The form factors G v {q 2 ), G M (q 2 ), G A {q 2 ), G P (q 2 ) are 
real functions of the squared momentum q 2 . They are 
parameterized as follows: 



Gv(q 2 ) 
G A (q 2 ) 



fJV 



1- 

Ml 
9A 



G M {q 2 ) 



9m 



1 - 



('-*) 



(11) 



The two form-factor cut-offs My and Ma are in general 
different and their values are expected to be of the order 
of 1 GeV like it is in the case of nucleon form-factors. 
As it will be discussed later the q 2 -dependence of these 
form-factors is not crucial for tritium /3-decay. 

The conserved vector current hypothesis (CVC) im- 
plies gv — 1-0. gM = —6.106 is calculated from the 
values of magnetic moments of 3 H and 3 He using the 
CVC hypothesis as well [27j . The axial coupling con- 
stant g A can be determined from the measured half-life 
of 3 H. The induced pseudoscalar coupling is given by 
the partially conserved axial-vector current hypothesis 
(PCAC) 



g P (q 2 ) = 2M t g A (q 2 )/(ml 



(12) 



TOtt is the mass of pion. 



For the spin-summed, Lorentz-invariant squared am- 
plitude we get 

~ \M\ 2 = l6(G F V ud ) 2 

spins 

x [G 2 v Vvv + GaGvVav + G\V AA + 



-G A GpP A p + GpPpp 



. r , n Vvm . r , n Vam 



2Mi 



2Mi 



n 2 'PMM 



T M 



4Mf 



with 



Vw = PefPui + PeiPuf ~ M t MfP ev , 



(13) 



(14) 



VAA = PefP„ l + + MiM f P ev , (15) 



Vav = 2 [P ef P vi - P ei P vf ) , 



(16) 



Pap = M f (miP„ t +mtP ei )-M i (miP„ f +mi,P ef ), (17) 



Vpp = -(P if ~M l M f )(P eu {m 2 e +ml) + 2mlm 2 e ) , 



(18) 



Mi [Peu{Pf " Mf) + P e f{P vl - 2P vf ) + P a P vf ] 

+M f [P e APf ~ Mf) + P ei (P vf - 2P vi ) + P ef P ui ] ,(19) 
Vam = 2{Mi + M f )(P e fP vi - P ei P vf ), (20) 



Vmm 



~^Pif (Pev(m 2 e + ml) + 2m 2 e ml) - M l M f m 2 e m 2 
+2P ei P ef (P ev + ml) + 2P vl P v f(P ev + ml) 

--MiM f P eu (3m 2 e + 3m 2 u + 4P el/ ) . (21) 

Here, Pu = (Pk ■ Pi) with k,l = i, /, e and v denotes 
the scalar product of two four-momenta. 

By neglecting the contribution from higher order cur- 
rents (terms proportional to Gm,p) we find 



1 



^ 'M| 2 = 16{G F V ud ) 2 



2 ^ 

spins 

x [(Gv + G A ) 2 (P e -Pf)(P,-P l ) 
+ (G v -G A ) 2 (Pe-P)(PvP f ) 
{-G 2 v + G A )MiMf{P e -P v )] . 



(22) 



The advantage of the presented formalism is that the 
squared Lorentz invariant amplitude is calculated ex- 
actly unlike in Ref. (2p| . where an assumption about 



its dominant constituent was considered. We note that 
for Gy = Ga = 1 the squared amplitude is proportional 
to (P e ■ Pf){P v ■ Pi), i.e., the structure is similar as, e.g., 
in the case of the muon decay. 

For the tritium /3-decay at rest the differential decay 
rate is 



dT = 



k F(z ' E ' } (5 £ H 

\ spins J 



2M % 

(2tt) 4 
(2tt) 9 



S^(P-P f -P e -P v ) 



d 3 p e d 3 p v d 3 p f 



2E P 2E U IE 



f 



(23) 




0.8 -0.6 -0.4 -0.2 



-y-m y [eV] 



The factor 1/2 in front of the squared amplitude stands 
for the average over the spin of the initial state. 

The subject of interest is the energy distribution of 
the electron. Hence, the integration over antincutrino 
and final nucleus momenta have to be performed in 
It requires calculation of the following integrals: 



K. = 
(A,,/)' = 



d 3 p f d 3 p v 
E f E v 



SP>(Q-P f -P v ) 



J^l^S^(Q-P f -P v )(P vJ y 



(24) 



-(25) 



d 3 p f d 3 p u 
E f E v 



^\Q-p s -p v ){p h y{p l r 

(26) 



with Q = Pi — P e and k, I — v,f. The details of integra- 
tions with results are given in the Appendix. 

The differential decay rate is found to be of the form 



dT 

dEL. 



-^(G F V ud ) 2 F(Z,E e ) Pe 
2n } 



M i I , „ M f 

{m 12 ) 2 y v M i 

x [gyTZvv + 9A9vTi-AV + g 2 A TZ A A+ 

+ gAgpTi-AP + gpTlpp 
+9v9mTZ-vm + 9a9mTZ-am + 9\i^mm\ , 

(27) 

where (m 12 ) 2 = M 2 +m 2 e ~2M l E e and y = E™ ax -E e . In 
the calculation we neglected q 2 dependence of the form- 
factors as for the /3-decay of 3 H the value of q 2 is rather 
small. Their consideration would lead only to small cor- 
rection factors, which are not sensitive to neutrino mass. 
We find not usefull to present here the explicit form of all 
Tli {I = VV, VA, AA, AP, PP, VM, AM, MM) fac- 
tors. Instead of that we conclude about their structure 
and importance. 

Our analysis showed that each term of IZi is propor- 
tional to (y + m v (M f + m v )/Mi) or (y + m v MfjMA. 
So, a common (y + m„Mf /Mi) can be put in front of the 
bracket in (f2"T)l by neglecting a small term m v /Mi. The 



FIG. 1: Endpoints of the relativistic Kurie plot [see Eqs. (|31[) 
and (|32[) ] of the tritium beta decay for various values of the 
neutrino mass: m v — 0, 0.2, 0.4, 0.6, and 0.8 eV. 



importance of different IZi contributions can be studied 
in the limit Mi = Mf, E e = m e and by making Taylor 
expansion in in m v , m e (m„ <C m e <C Mi). The leading 
terms of different IZi (without the common factor) are 
as follows: 



VV : m e Mi, AA : 3m e Mi, AV : 2ml, 

AM : 2m 2 e , 
4 



VM: l< 
2M, 



MM: ±< 
16 Mf 



Tfl \ HI 

AP: 2m e M^, PP : -m^—^r 
mi 2 Mfmi 



(28) 



From their comparison we conclude that the contribu- 
tions coming from higher order terms of hadron current 
to the decay rate of 3 H can be neglected. 
Then we have 



dT 

dEL. 



1 

2^3 



(G F V ud ) 2 F{Z,E e ) Pe 



Mf 



'y y + 2m v 



M f 



(mi2) 2 ^V ' " Mi 

, 2 / MA M 2 {E 2 -m 2 e ) 

{gv + 9A) 2 y[y + m v 7 > 

{gv + 9A) 2 (y + m v 



X (y + M f 



Mi J 3(m 12 ) 4 
Mf +m v {M l E e -m 2 e ) 

Mi m{ 2 
M f + m v .{Mf-MiE e ) 



Mi 



(9v ~ 9a) Mf y + 



(M f + M v ) 



M, 



(M t E e - ml) 



(mi 2 ) 2 

+ (gv ~ 9a) 2 E £ [y + m v 



Mf 
Mi 



(29) 



The first term in the brackets in (f2"9"| . which is quadratic 



■5 



in y, plays a sublcading role. By keeping only the dom- 
inant contributions and by introducing a mass scale pa- 
rameter M instead of the Mj and Mf, we get 



7T"1 -J 

— ~ ^(G F V ud ) 2 F(Z,E e ) Pe E e (g 2 v + 3g 2 A ) 



x y/y{y + 2m v ) (y + m v ) ■ 



(30) 



For the relativistic form of the Kurie function we can 
write 

K(y) = B (y/i (y + 2m v ) (y + m„)j V ' (31) 



with 



GpV u< i , 

B = . = \hy 



(32) 



The unknown coupling constant q_A of the hadron current 
is fixed to the half- life of 3 H [22, [23[ with result qa = 
1.247. This value coincides well with that of the axial- 
vector coupling of the nucleon (see previous section). We 
have B = 3.43 x 10~ 6 GeV~ 2 . 

By comparing the Kurie function in (|3ip and (|32p with 
the commonly used one ^ we find that they are equal 
if y is replaced with (Eq — E e — to„) and |Mgt| 2 = 3 
is assumed. This confirms what was generally expected, 
namely that the relativistic effects are small corrections 
to the results known in the traditional method due to 
a small Q-value of the /3-decay of tritium. However, it 
was not clear yet whether the recoil of the nucleus, which 
value is 3.4 eV for maximal electron energy, affects the 
cndpoint spectra, if sub eV mass of neutrino is measured. 
Within the considered EPT of /3-decay of tritium we find 
that there is no significant modification of the shape of 
the electron spectra close to the endpoint due to the nu- 
clear recoil. 

In Fig. [1] we show a relativistic Kurie plot for the j3- 
decay of 3 H versus y — E max — E e near the endpoint. 
Special attention is given to the effect of a small neutrino 
mass (m„ = 0.2, 0.4, 0.6 and 0.8 eV). We see that the 
Kurie plot is linear near the endpoint for zero neutrino 
mass (m u — 0). However, the linearity of the Kurie plot 
is lost if the neutrino has a non-zero mass. Deviation 
from a straight line depends on the magnitude of neu- 
trino mass m v . Though, there is no difference with the 
previously known dependences, it is worth to stress that 
in this case the relativistic form of the /3-decay Kurie 
plot is used, which also takes the nuclear recoil (~ 3.4 
eV) into account. 



IV. CONCLUSION 

The neutrino absolute mass scale, which is very impor- 
tant for particle physics as well as for cosmology and as- 
trophysics, cannot be resolved by oscillation experiments. 
A way of the direct determination of the neutrino mass 



scale in laboratory experiment is the investigation of the 
kinematics of tritium /3-decay. 

The KATRIN experiment [H, [H, E3] , which is under 
construction, will be able to reach a sensitivity of neu- 
trino mass in the sub-eV range. In connection with that 
there is a request for a highly accurate theoretical de- 
scription of the electron energy spectrum. 

In this paper we derived the relativistic form for the 
/3-decay endpoint spectrum in the elementary particle 
treatment of weak interaction. The considered formal- 
ism follows from the analogy between 3 H ( 3 He) and the 
neutron (proton) having the same spin-isospin proper- 
ties. It allowed us unlike in Ref. [2C| to determine the 
squared /3-decay amplitude more accurately. In addition, 
we found that the higher order terms of the hadron cur- 
rent can be neglected without affecting the dependence 
of the Kurie plot on the electron energy and the neutrino 
mass. By comparing the relativistic and previously used 
Kurie functions a good agreement between them was es- 
tablished. 

We acknowledge the support of the EU ILIAS project 
under the contract RII3-CT-2004-506222, the Deutsche 
Forschungsgemeinschaft (436 SLK 17/298) and of the 
VEGA Grant agency of the Slovak Republic under the 
contract No. 1/0249/03. 



APPENDIX A 

Here we outline the calculation of integrals over neu- 
trino and final nuclear momenta. 

Integration of JC: 

The integration is performed by choosing Q = (Q n ,0), 
i.e., the rest frame connected with the center of mass of 
antineutrino and final nucleus. We have 



K. 



J -jjr S(Q -Ef- E v ) PfdE f dn f 



(Al) 



with E f = {ml - Mj + Ej) 1 / 2 . By using 5(f(x)) 
S(x - x Q )/|/'(xo)| we find 



K = 2ir 



^[Ql - (M f + m„)2][Qg - (M f 

Ql 



(A2) 



We replace Q 2 with Q 2 and write K. in the Lorentz in- 
variant form 



K = 2?r 



V[Q 2 - (Mf + m„)2][Q2 - (Mf - m v f 



Q 2 



y(y + 2m v ^f- 



(m 12 ) 2 



(A3) 
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Integration of {C v ) p : 
The integral 



(A4) 



is a second rank tensor 



(M,/)" = Cg pa + DQ"Q° 



(A7) 



(A8) 



can be written as 



(£„)' = AQP. 



(A5) 



Here, A = A(Q 2 ) is a scalar function of Q 2 . By multiply- 
ing (C u ) p with <2 p the constant A(Q 2 ) can be determined. 
Then we get 



(A6) 



Integration of (AC/) P<T : 
The integral 



CM,,)" = 



E f E v 



5W{Q-P f -P v ){P v )»{P f )° 



where C = C(Q 2 ) and £) = D{Q 2 ) are scalar functions 
ofQ 2 . 

By multiplying {MvfY" with g^ v and with Q p Q a a set 
of two equations is formed. By solving them we find 



CM,/) 



/ItT 



{Q-P V ){Q-Pf)\ 1 



(Q ■ P v ) (Q ■ Pf)\ Q P Q° 



3Q 2 ' 
(A9) 



The remaining integrals {C v ) p , {M vv ) pa , (A/>/) p<T can 
be calculated following the scheme given above. 
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